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Abstract 

In a recent publication (Bertulani et al., PLB 650 (2007) 233 [1]) the validity of 
analysis methods used for intermediate-energy Coulomb excitation experiments was 
called into question. Applying a refined theory large corrections of results in the 
literature seemed needed. We show that this is not the case and that the large 
deviations observed in [1] are due to the use of the wrong experimental parameters 
in that publication. We furthermore show that an approximate expression derived 
in [1] is in fact equivalent to the theory of Winther and Alder, an analysis method 
often used in the literature. 
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1 Introduction 



Coulomb excitation (CE) of radioactive heavy-ion beams on high-Z targets at 
intermediate to relativistic beam energies (E/A > 30 MeV) is an important 
and now widely used method to study the structure of nuclei furthest from 
stability and a wealth of experimental information was obtained in laboratories 
worldwide [2,3,4,5,6] since the method was first applied in 1995 [2]. Even 
at these high beam energies it is possible to suppress the influence of the 
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strong interaction in the excitation process, by rejecting events with small 
impact parameter, and to extract largely model independent nuclear structure 
information. 

Experimentally, usually the de-excitation 7 ray yields after CE are measured 
and from these the excitation cross sections are deduced. It is important to 
note that the deflection angle of the outgoing scattered particle is restricted 
either by the experimental setup |2,3,7,8j or during the off-line analysis of 
the data, which are recorded event-by-event |9j. The so imposed maximum 
scattering angle 9 max corresponds to a certain minimum impact parameter 
b m i n which is chosen to exceed the sum of the radii of projectile and target by 
several femtometer, ensuring the dominance of the Coulomb interaction in the 
excitation process. The angle integrated CE cross section depends of course 
strongly on the chosen maximum scattering angle 9 maxi which is therefore 
an important parameter in the analysis. For further details review articles 
[T0~jll|ll2j or the original references should be consulted. 

In order to extract the relevant quantities, namely the reduced transition 
probabilities B(ir\) and in particular £>(E2; 0+ — > 2+) values, various analysis 
techniques are used: 

• distorted-wave theory (using the ECIS code [13]) used e.g. in [2|14II15II16] . 

• the semi-classical theory of Winther and Alder [171 used e.g. in [31I91T8] . and 

• the virtual photon method [T9] used e.g. in [2PjJ_j 

More than one method has been applied to certain nuclei and consistent results 
were obtained, see e.g. [2|14|21|22j . Recently, a compilation of data obtained on 
stable nuclides was published and a comparison of the -B(E2; 0+ — > 2±) values 
obtained by various methods was given [23J. It was concluded that high-energy 
CE is at least as precise and accurate as so-called standard methods such as 
direct lifetime measurements, low-energy "safe" CE, and electron scattering. 

Recently, however, a new approach was proposed [21] and previously published 
results were called into question in [TJ claiming that either only non-relativistic, 
or only partially relativistic, expressions are used or the effect of the orbit of 
the nuclei, due to their Coulomb repulsion, is not fully accounted for. Using 
their new model, changes in the calculated cross sections of up to 30% were 
reported with respect to results from the literature. 

We would like to dispute that assertion and show that the previously published 
results 13,25,91 are indeed correct and consistent. 



The virtual photon method is in fact identical to the semi-classical theory of 
Winther and Alder. Only here the excitation cross section is represented as a product 
of a virtual photon number and a (virtual) photo absorption cross section. 
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2 Cross sections 



In pp Bertulani et al. studied the CE cross sections at intermediate energies 
(10-500 MeV/w) using their model, published earlier in [24], which aims to 
correctly describe the CE process at all beam energies. In addition, since the 
orbital integrals appearing in the expressions for the excitation amplitudes 
have to be determined numerically, an approximate expression for the CE 
cross section was derived, called a app or a app in pQ, which is a simple expres- 
sion containing modified Bessel functions. Using these expressions and using 
the numerically determined orbital integrals it was found that the calculated 
excitation cross sections for a given excitation strength are smaller by 10-30% 
in comparison to values from the literature [3.9.25J and consequently that the 
published £?(El)f and _B(E2)| values are in fact too small by 10-30%. 

Since this would be a substantial alteration we tried to reproduce (i) our 
original results in [3|9|25] and (ii) the results of Ref. pQ using Eq. (5) of pQ. 
Our findings are listed in table [1] and are displayed in figure [fl The number 
shown in column 1 of the table and plotted on the abscissa in the figure 
corresponds to the "Data set" as given in pQ. Contrary to Bertulani et al. we 
essentially found exact agreement between the published cross sections (circles 
in the figure) and the results of our own calculation using Eq. (5) of [1] (filled 
diamond in the figure). However, the results reported in pQ (triangles) deviate 
significantly from those. 

In order to understand the differences we performed a calculation using the 
beam energies and scattering angles listed in table 1 of pQ. We would like 
to point out that at least for data sets number 7-12 (original publications 
]3 9j) the given values correspond to the incident beam energy and maximum 
laboratory scattering angle. If we assume, however, that the listed beam energy 
is the average beam energy in the mid-plane of the CE target^] and the angle 
is the scattering angle in the center-of-mass system, we can in all cases, except 
for 11 Be, reproduce the values given in pQ (see open diamonds in figured]). 

We therefore conclude that mentioned discrepancies are due to an interchange 



2 To determine the effective (theoretical) cross section a the beam energy dependent 
cross section a(Eb) should be averaged over the target thickness dt 



where Ef ) (x) is the beam energy at target depth x. Due to the almost linear depen- 
dence of the cross section on beam energy within the small energy loss in the target 
it is sufficient to use the beam energy in the center of the target in the calculation. 
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Fig. 1. Experimental CE cross sections as listed in table [TJ The number on the 
abscissa corresponds to the "Data set" of table [TJ and is identical to the one given 
in pp. The circles show the experimental cross sections a exp and the filled diamonds 
the cross sections calculated by us according to Eq. (5) of pp (cr' app in rows with "*" 
in table [1]). The open diamond shows the result of our own application of Eq. (5) 
of pp using the input parameters from the rows without a "*" , i.e. the ones used in 
PP. The triangle shows the result published in pp. 

of center-of-mass and laboratory angles and the use of the incident beam 
energy instead of the mid-target beam energy. 

We furthermore found that the expressions given for cr app in pQ are in fact 
analytically identical to the one of Winther and Alder [T7] (see Appendix). 
Since for all experimental cross sections listed in table 1 and shown in figure 
1 the semi-classical theory of Winther and Alder was employed during the 
analysis of the data, the a app results must be identical to the original results 
a exp provided the proper input parameters are used. 

This conclusion is in fact supported by an earlier work of Bertulani et al. [23] 
where it is shown in figure 3 (curve RR) that for the case of an E2 excitation 
the cross section determined within the Winther and Alder theory deviates by 
less than 5% from the exact treatment for beam energies above 30 MeV/u. 



We would also like to note that the analysis of CE data using a simple max- 
imum scattering angle is not applicable for all experiments in the literature. 
For instance for the analysis of the CE of 32 Mg [2] the use of a sharp cutoff 
6 max is meaningless ("Data set" 6 in pp). Here the particle detection efficiency 
depended (smoothly) on the scattering angle with maximum values of about 
80%. In this case the differential cross section must first be multiplied by the 
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Table 1 



Coulomb excitation cross sections are shown in the three columns on the right-hand 
side with (J exp being the experimental CE cross section, a app the cross section from 
reference pQ and a' app the cross section calculated by us using Eq. (5) of pQ. The 
B^EX)^ values are listed in units of e 2 fm 2 \ Values in bold were taken verbatim from 
PQ and the "Data set" number # corresponds to the one used in pQ. In the rows 
starting with "*" the correct center-of-mass scattering angles 0'f^ ax and mid-target 
(m-t) b earn energies B^ eam are listed; an agreement of (Jexp and cr ! a pp is evident for 
all cases. The rows not starting with "*" list the values as given in pQ and a' app is 
calculated by us using these values. In most and a' app agree, but a eX p and 

a apv do not. 

experimental particle detection efficiency before integrating over the scatter- 
ing angle. This can only be done when the efficiency curve is known, which is 
shown in figure 2 (lower panel) in [2j. 

We furthermore found the following errors in pQ: 

• There is a misprint in Eq. (5). For the E2 excitation cross section the power 
of £ 7 should be 2 (E^) and not 3 (£?). 

• The wrong citation was used for 38 > 40 > 42 S and 44 ' 46 Ar. It should be Scheit et 
al. [3] and not Chromik et al. [20] . 

• An Ml excitation was considered for 17 Ne, but the excitation process pro- 
ceeds through an E2 transition with only a very weak Ml component. The 
decay, though, is an almost pure Ml transition. 
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• The deviation of the cross sections published in [TJ (triangle in figure [T]) for 
n Be ("Data set" 3) seems to be due to a forgotten factor of e 2 = ahc ~ 
197/137^ 1.44 in jlj. 



3 Conclusion 

We conclude that the theory presented in [TJ24] does not represent a significant 
improvement over analysis methods used so far. In fact, a good approximation 
[TJ to the full theory is equivalent to the theory of Winther and Alder |17j . 
which is widely used in the analysis of CE data at high beam energies. We 
furthermore conclude that the large differences of up to 30% reported in [TJ 
are likely solely due to an interchange of center-of-mass and laboratory angles 
and the use of incident beam energy instead of the mid-target beam energy in 
the calculation of the CE cross sections. 

We nevertheless agree with the authors of [TJ that a refined theory which 
includes relativistic effects and properly treats the orbit of the projectile is 
needed and should be applied in the future, when high-precision results be- 
come available at next-generation radioactive nuclear beam facilities. The cur- 
rently available data will be hardly affected due to rather large statistical and 
systematical uncertainties. 
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A Coulomb excitation cross sections 

In the following o B denotes the expression for a app from Bertulani et al. [TJ 
Eq. 5] and aw /a the result given by Winther and Alder [17] . For a given 
multipolariy A these two expressions for the CE cross section differ explicitly 
only in the expressions containing the modified Bessel functions K n and we 
will consider the ratio of the two expressions leaving out explicit common 
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pre-factors. In particular for the E2 case: 



(E2) 



a 



W/A 



E \G E 2»(f3- 1 )\g»(0 

fj,=-2,...2 



(A.1) 



where the functions and GVa^ are given in pTTl Eqs. 2.16 and 3.4]. The former 
contain the Bessel functions K^+x and the latter is a polynomial function of 
f3^ 1 . After a lengthy, but straight-forward, calculation one obtains: 



+ 3 7 - i )K 2 2 - (8 7 ~ 2 + W)K X K 2 + £7 #3 - 4 7 ~ 2 i^: 



(A.2) 



Using the recurrence relations for the Bessel functions K n [26] 



K — K 2 = and 



-K x + K 3 = -K 2 (A.3) 



the Bessel functions K„ with /i > 1 can be expressed in terms of and K\ 



K 2 = -K, + K 

4 /8 \ 4 

K 3 = ^2 + K, = f - + lj KTx + -K . 



(A.4) 
(A.5) 



The result is 

£|gW(0~ 1 )%m(0 



4tt 2 
75^ 



/3 4 e 2 ^o 2 + + 4 7 ~ 2 )^ 2 + 2£(1 + 7" 2 ) 2 K K 

(A.6) 



i > 



where /3 4 + 47 2 = (1+7 2 ) 2 was used. The last result (IA.6[) is identical to 
the numerator of (1A.1I) and therefore 



(£2) _ (E2) 



(A.7) 



Similarly for the El case: 
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(El) 32 
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(El) 
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(A.8) 



7 



and 

E | W/nlVo = ^(V 2 - iwi - kd + a-\-K\ + k\ - H^)) 

(A.9) 

Using the above recurrence relations for the Bessel functions we obtain 

E IWr 1 )! 2 ^) = ^^(^o 2 - Kl + -^K oKl ) . (A.10) 
The last result (lA.lOj) is identical to the numerator of (IA.8I) and therefore 



A E1) = ^f A (ATI) 



also for the El case. 
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